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Computer modeling is reported about the creation and pinning of a magnetic vortex-antivortex (V- 
AV) pair in a superconducting thin film, due to the magnetic field of a vertical magnetic dipole above 
the film, and two antidot pins inside the film. For film thickness = 0.1£, k = 2, and no pins, we find 
the film carries two V-AV pairs at steady state in the imposed flux range 2.10<f>o < <f> + < 3.0<I?o, 
and no pairs below. With two antidot pins suitably introduced into the film, a single V-AV pair 
can be stable in the film for $ + > 1.35>o- At pin separation > 17£, we find the V-AV pair remains 
pinned after the dipole field is removed, and, so can represent a 1 for a nonvolatile memory. 

PACS numbers: 79.60.Bm, 73.20.Dx, 74.72.-h 



X 



I. INTRODUCTION 

The electronic properties of superconducting devices 
are critically influenced by the motion of vortices (Vs) 
within the superconductor (SC). Pinning of Vs is respon- 
sible for critical current and magnetic hysterisis in the be- 
havior of a SC. 1 Any defects, such as magnetic dots, twin 
boundaries 2 , local spots with lower critical temperature 3 
T c or variable thickness, 4 can act as pinning centers. 

Ferromagnetic particles, fabricated onto a supercon- 
ductor, have many effects of fundamental interest. Van 
Bael et al. 5 considered a ferromagnetic dot array with 
magnetic dipole moments parallel to the surface of su- 
perconducting substrate. They showed that a flux lattice 
(FL) is pinned at the opposite poles of the ferromagnetic 
dots (FDs), where flux quanta of opposite signs are in- 
duced by the stray field. Nozaki et al. 6 investigated the 
effect of a spatially modulated local field by a ferromag- 
netic dot array on a Nb film. Since the magnetic field line 
from the ferromagnetic dots supresses superconductivity 
where it penetrates, these magnetic dipoles (i.e. FDs) 
can effectively create vortices and antivortices, and also 
have a pinning effect on those vortices. 

Milosevic and Peeters 7 solved the Ginzburg-Landau 
(GL) equations for the case of a ferromagnetic disk (FDk) 
on top of a thin superconducting film. They found that, 
as magnetic field lines penetrated a localized region of 
the film and then returned to the outside, a giant an- 
tivortex (AV) formed under the FDk, and several Vs lo- 
cated themselves in a ring outside the region in the film 
below the FDk. Priour and Fertig 8 performed a sim- 
ilar Ginzburg-Landau study for the vortex states in a 
superconducting thin film, subject to the magnetic field 
of an FD array, with the dipole moments oriented per- 
pendicular to the film. By varying the dipole moments, 



they changed the number density of V-AV pairs, with the 
system going through various vortex configurations that 
even break dot lattice symmetries. 

Submicron-sized holes, or antidots, are well-known as 
strong pins. Columnar defects created by the heavy-ion 
irradation technique are known as efficient pinning cen- 
ters, strongly enhancing the critical current density. 9 ft 
has been shown 10 ' 11 that a small antidot induces nega- 
tive potential energy around its center, attracting a V 
as a singly-quantized flux line. Budzin 12 has shown that 
for large enough antidots, multiple flux quanta can be 
trapped. 

Takezawa and Fukushima 13 presented a simulation for 
the pinning force on a V moving around a square an- 
tidot. They solved the GL equations, and the dot size 
was varied to find the maximum pinning force. Priour 
and Fertig 14 also performed a numerical study of the be- 
havior of Vs in the presence of an array of antidots, in a 
thin superconducting film, using the GL equations. They 
showed that when a V moves within a critical distance, 
d c , of an antidot, the V core deforms and extends to the 
antidot boundary, and the associated supercurrent and 
magnetic flux spread out and engulf the antidot. 

Pearl 15 solved the London equation for vortices moving 
in a thin film, and showed that the vortices have a longer 
range for their interaction force than in the bulk SC. This 
leads to the concept of an effective penetration depth 
X e ff = A 2 /2d, which is much longer than the magnetic 
penetration depth A in thin films where d <C A. 

In this paper we report computer modeling of: (1) the 
creation of 2 and 3 V-AV pairs, but not one pair, in a 
thin uniform superconducting film with a vertical mag- 
netic dipole above it. The dipole is made of two magnetic 
monopoles of strength ±g m , separated by a distance di, 
with the positive monopole closer to the film at a height 
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z m above the film; (2) the creation and pinning of a single 
V-AV pair in the film with two antidot pins introduced 
in the film at appropriate locations and separation; and, 
(3) the required V-AV separation to keep the pair pinned 
at two antidots even after the dipole field is removed. 16 
The motivation for this investigation is to investigate the 
scientific feasibility of a non-volatile memory by using a 
magnetic dipole to create a V-AV pair in a supercon- 
ducting thin film, that relies on two properly positioned 
pinning centers in the film that pin the V-AV pair and 
prevent them from annihilating each other after removal 
of the dipole. However, we do not attempt to analyze 
whether such a memory device can compete with other 
existing memory devices in cost and speed, we leave such 
concerns to future studies. 



Using a link variable approach 20 ' 21 on a staggered grid 
over f2 shown in Fig. 1, the semi-discrete equation is: 
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where d = d(x,y) denotes the film thickness and 
(h x ,h y ) the spatial increments in (x, y)-directions. Note 
that the subscript of each variable (except x and y in h x 
and h y ) denotes the location at which it is defined on the 
grid (see Fig. 1). 



II. THE SIMPLIFIED GINZBURG-LANDAU 
MODEL 

For a thin superconducting film with d < £ (the coher- 
ence length) and A, we neglect variations of all quantities 
in the thickness direction of the film. 17 After nondimen- 
sionalization, the GL equations take the form: 
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FIG. 1. The staggered grid arrangement for cell nodes 
P,E,W,N,S and faces e,w,n,s. 
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Distances are measured in units of £; the vector potential 
A is in <I>o/27r£, where $o = hc/2e is the flux quantum; 
the magnetic field is in H C 2 = $o/2< 2 = V2kH c , where 
k = A/£ is the Ginzburg-Landau parameter. The London 
gauge V • A = has been employed for A. The dimen- 
sionlcss gauge-invariant free energy functional is then: 
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The vector potential A is obtained by a fast Fourier 
transform (FFT) technique. 22 The 2-dimensional super- 
current j 2 D = (jx,j y ) on the grid is given by: 



ixe = — Pp6e - ep* E )cos(A e /i x ) - (*p$ E + e P e E ) sin (A e hx)] , 
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where $ and <d are the real and imaginary parts of ^, 
i.e. * = $ + iQ. 

In this study, we have assumed k = 2, h x = h y = 
0.25, and At = 0.05 and 0.1 for spatial and temporal 
discretization. For more details sec 23 . 



where the ^/-dependent terms are integrated inside the 
film only. The first GL equation is solved by a relax- 
ation method 18 for a pseudo-time-dependent equation 
= — jpr), in a square sample with periodic bound- 
ary conditions 19 . We are primarily interested in finding 
steady states (note that the London gauge V • A = 
is a good choice only at steady state). The relaxation 
method provides an effective means to solve the nonlin- 
ear PDE, where the order parameter changes in a pseudo- 
time, while the vector potential and the supercurrent are 
updated simultaneously. 



III. MAGNETIC PENETRATION IN A FILM OF 
CONSTANT THICKNESS 

The film is subject to a non-uniform magnetic field 
from a magnetic dipole above the film, and has pinning 
centers such as holes (antidots) of various sizes and lo- 
cations. Since dfu m = 0.1£ <C A, the film is virtually 
transparent to magnetic penetration in the thickness di- 
rection, and is "soaked" in the applied magnetic field 
(Fig. 2(a)). 
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FIG. 2. (Color online) (a) A superconducting film under 
a magnetic dipole. Stray field lines penetrating the film cre- 
ate vortex-antivortex pair(s). They may be pinned at pin- 
ning centers such as holes (antidots) made in the film in a 
proper configuration. Since dfu m — 0.1£, the film is virtually 
transparent to magnetic penetration in the thickness direc- 
tion, and is "soaked" in the applied magnetic field. In the 
film, the induced magnetic fields by the supercurrents incur 
only negligible change in the applied magnetic field, (b) The 
3-dimensional plot of $ + /g m due to a single dipole. <E> + /q m 
decreases as z m increases, and increases as d; increases. 



Stray field from the magnetic dipole passing through 
the film creates V-AV pairs. Since the Vs and AVs 
are topological excitations of opposite quantum numbers, 
there must be an equal production of them in the film. 
Field lines change polarity in the film so that the total 
flux penetrating the film is zero. We define the integral of 
the downward part of the dipole field in the plane of the 
film, i.e., <& + = fxx ^ n Hi/ • nda , where n is its down- 
ward normal, as the imposed magnetic flux in the film by 
the magnetic dipole (see Fig. 2(a)). $ + = $ + /$o pro- 
vides a rough guide on how many V-AV pairs the dipole 
attempts to create in the film. 

For a single dipole, the magnetic flux can be written 
exactly as: 

$+ = -i r 2npH z (p,0)dp 



2tt 



(7) 



where we have used the following expression for the z- 
component of magnetic field: 
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In the above p = x + y on the film, z m = z m + \df, 
and p c is the radial location on the film where the z- 
component of the magnetic field changes its polarity. It 
is given by the exact expression p\ — Zm 3 (z m + d/) 4 / 3 + 

Zm(z m + di) 2 / 3 , when due to a single dipole. But for a 
periodic array of dipoles, we have to evaluate it numer- 
ically. (We make the thin film approximation by com- 
puting <E> + at the central layer of the film, which is at 
z = 0.) As p c increases the magnetic field lines reach 
farther outside. However, to localize the field lines in- 
side the computational domain, we cannot use z m and di 
that are too large. Instead, we have to consider z m and 
di that are much smaller than the domain size. For this 
reason we have not systematically explored their individ- 
ual influence on the results obtained. Instead, we have 
explored the effect of the most important parameter $ + 
which can be achieved by varying q m , z m , and/or d/. 

In addition, to consider a finite domain with periodic 
boundary conditions, we have to use a periodic array the 
dipoles (49 x 49 dipoles were superposed to simulate the 
periodic array). The magnetic flux <f> + calculated for 
the periodic array of dipoles differs from that using a 
single dipole by less than 5%, in the range of z m and di 
considered. 

Figure 3 shows magnetic penetration into a thin film 
of thickness df — 0.1, with no antidot pins. The compu- 
tational domain is 16£ x 16£, with periodic extensions of 
this geometry. We fix z m = 0.8£, q m = 20<I>o, an d vary 
di to change $ + . At $ + = 3.0 we find three V-AV pairs 
[Fig. 3(c)], and for 2.10 < $+ < 3.0 the film has two 
V-AV pairs [see Fig. 3(b), for $+ = 2.50]. 
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The AVs accumulate near the center of the him, (i.e., 
directly below the dipole), to form a giant AV, whereas 
the Vs spread outward symmetrically because they repel 
each other and, although they are attracted to the giant 
AV, they prefer to exist in the region where the magnetic 
held points upward. 

Plots of the phase of the order parameter give a clearer 
picture. Thus, in Figs. 3 (a), (b) and (c) respectively, 
the density plot is at the top and, for (b) and (c), the 
corresponding phase plot is at the bottom, and the phase 
varies by 2tt around the core of each (singly-quantized) V. 
It varies by 47r(67r) in the opposite direction around the 
core of a doubly(triply)-quantized giant AV. The phase 
plots show that the Vs and AVs are paired one to one. 

However, no V-AV pairs are formed for $ + < 2.10. In 
Fig. 3 (a) ($ + = 2.05), the magnetic held penetrates a 
large area in the him, just below the dipole. The area 
near the center of the him is largely penetrated due to 
the high value of the magnetic held around the center. At 
steady-state, no vortex forms and the screening supercur- 
rent (shown at the bottom of Fig. 3) is circulating only 
in one direction, as required by the theory of London 24 
for the Meissner state under a dipole held. Fig. 3(d) is a 
phase diagram showing the steady state Gibbs free energy 
V = V I A a ( k2A ' J 2D ~ 5l*l 4 ) dA - As * + increases, 
the rising magnetic penetration (vortices) also lifts G/V 
by increasing the electromagnetic interaction energy and 
lowering the condensation energy. The dashed lines sepa- 
rate regions of different number of V-AV pairs: 0, 2, and 
3 pairs (denoted by Roman numerals). Also, the free 
energy curve steepens going from 2 to 3 V-AV pairs. 

Milosevic and Peeters 7 suggested that the number of 
V-AV pairs increase with $ + in a non-integral increment 
of A$ + = 1.073 when they used an FDk dipole of radius 
r c i = 1£ and thickness dd = 1£. In their study, more 
V-AV pairs tend to be created as the radius of the FDk 
decreases for a given dipole strength. Our dipole, consist- 
ing of a magnetic monopole-antimonopole pair separated 
by a small distance, corresponds to an extremely thin 
magnetic tip above the him and produces a highly fo- 
cused magnetic held. Thus, our A$ + = 0.85 between 
Fig. 3(a) and (c) agrees qualitatively with their con- 
clusion, implying that it becomes easier to create V-AV 
pairs as the lateral size of the dipole decreases. Fur- 
thermore, in Ref. 7 the (dimensionless) dipole moment 
m/7J C 2^ 3 ~ 11 — 15 is necessary to create 2 V-AV pairs for 
Rd = 3.5£, whereas our dipole moment q m di ~ 5.2 — 8.3 
with $ + covering about the same percentage of the him 
area, creates the same number of V-AV pairs. Thus, our 
point dipole (magnetic tip) is about twice as effective at 
creating 2 V-AV pairs than the FDk in Ref. 7. This we 
believe is because our dipole held is highly non-uniformly 
distributed within the area of radius p c . 

The magnetic dipole considered here generates a highly 
inhomogeneous magnetic held in the superconducting 
him. The very strong downward held, in a small region 
just below the dipole, holds all AVs to form a giant AV. 




(a) (b) (c| 
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FIG. 3. (Color online) The steady-states of magnetic penetra- 
tion into a uniform film with no pins under the external field from a 
magnetic dipole. The top is the contour plots of the superelectron 
density (blue = 0, red = 1), and the bottom is the corresponding 
phase plots (blue = 0, red = 2n) of the order parameter except 
the vector plot for supercurrent in (a). Here <&+ = 2.05 for (a), 
l>+ = 2.50 for (b), and i+ = 3.0 for (c). The film carries three 
V-AV pairs in the steady state shown in (c), two pairs in (b), and 
no pairs in (a), (d) shows the phase diagram of the steady state 
Gibbs free energy £ = ± J A (k 2 A • j 2 D ~ |l*| 4 ) dA. As l>+ 
increases, the rising magnetic penetration (vortices) also lifts G/V 
by increasing the electromagnetic interaction energy and lowering 
the condensation energy. The dashed lines separate regions of dif- 
ferent number of V-AV pairs: 0, 2, and 3 pairs (denoted by Roman 
numerals). Also, the free energy curve steepens going from 2 to 3 
V-AV pairs. 
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The much weaker upward field outside this region is not 
effective at pulling the V's away from the giant AV. A mu- 
tual repulsion between the V's also plays an important 
role in determining the final V arrangement. 7 Indeed, the 
absence of such a repulsion in the case of a single V-AV 
pair may account for the fact that we cannot obtain a 
single V-AV pair in any <£> + range. In Ref. 7, this lack 
of repulsive force was compensated by the effect of the 
finite-size ferromagnetic disk. 

Our computer modeling is based on a relaxation 
method with a pseudo-time. 18 Thus, computed transient 
states can provide some qualitative information about 
the way the system evolves in real time. In Fig. 4 two V- 
AV pairs are formed at $+ = 2.80. Initially, three V-AV 
pairs are created. Afterwards a pair of V-AV annihilate 
each other and two pairs remain in the steady state. This 
is more discernable in the phase plots (cf. Fig. 4.) 

We find that for 2.20 < <l + < 3.0, three V-AV pairs 
are created first. Afterwards, one V-AV pair annihilates 
itself, leaving only two pairs in the final steady state. The 
steady-state configuration shown in Fig. 3(b) belongs to 
such a situation. (This annihilation has not been ob- 
served at <f> + = 2.20, because the film carried two pairs 
throughout, suggesting it was close to the lower boundary 
of the two V-AV pair region, below which no V-AV pairs 
exist in the steady state.) Since our relaxation method 
begins with full penetration of the dipole field, a high 
peak field just below the dipole initially causes a larger 
number (3) of V-AV pairs to be created. Eventually, as 
the Vs separate from the AVs and move into a weaker 
upward field region, the system settles for only two V- 
AV pairs over the given $ + range. In addition, the vor- 
tices arrange themselves diagonally at the farthest pos- 
sible distance conforming to the given sample geometry 
and the periodic boundary conditions. 

The transient states discussed here are, of course, ob- 
servable at best in dynamic measurements only, since not 
only they are not thermodynamically stable states, we 
believe there are not even any metastable states involved 
at any intermediate stage. Some transient states can be 
quite long lived though, when they are nearly symmetric, 
such as the three V-AV pair state shown in Fig. 4 be- 
fore one V-AV pair eventually annihilate each other, with 
two remaining V-AV pair then rearrange their positions 
to produce the final stable state. 

IV. MAGNETIC PENETRATION IN A FILM 
WITH PINS 

When a pair of antidot pins are present in the film at 
suitable positions, we find that one V-AV pair can be 
formed as low as $ + = 1.3. Under this condition a pair 
of V-AV is created and then separated, each of which 
is attracted to an antidot by its pinning force. In the 
final state, an AV is pinned at an antidot placed near the 
center of the film (PI), and a V is pinned at an antidot 




t=0 t = 10,000 t= 19.000 




10 20 30 40 50 6 ID JO 30 40 50 60 ID JO 36 40 50 65 



t= 21,000 t =40,000 t= 70,000 



FIG. 4. (Color online) Pseudo-time sequence of magnetic 
penetration into the film with no pins under the external field 
from a magnetic dipole. At each time a pair of plots is shown, 
with the density plot at the top and the corresponding phase 
plot at the bottom. Here <3> + = 2.80 and is in the interval 
2.15 < <I> + < 3.0 discussed in Fig. 3. Three V-AV pairs 
initially appear in the film early (t=10,000), then one pair 
annihilates itself (t=19,000 and 21,000). This process is more 
discernible in the phase plot. Finally, the film carries two 
V-AV pairs at steady state. 
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FIG. 5. (Color online) For <§+ = 1.90 the steady-states of one 
V-AV pair creation in a film with two pins, under the external field 
from a magnetic dipole. In (a) and (b), two pins, PI and P2, are of 
the same size and placed the same distance apart (17/4)\/2£. The 
V-AV pair are nicely pinned in (b). However, a small deviation 
from the optimal pin position in (a) prevents full pinning (the AV 
stays near the edge of PI without being fully pulled in). Changing 
the pins in (a) to the bigger pins in (c) still does not achieve full 
pinning. 

some distance away along the diagonal (P2). 

The effects of dipole field and pinning force on the final 
location of the vortices is shown in Fig. 5. Here $ + = 
1.90, and for the two antidot pins we use square holes 
clipped at their corners so they have octagonal shape. 
The "diameter" d p of these approximately round pins is 
defined to be the edge length of the undipped square 
hole. In each case of Fig. 5 (a) and (b), two pins with 
dp = 2£ are placed at the same distance (17/4)\/2£ apart 
along a diagonal. 

For (a), the center of the pin PI is v2£ to the lower left 
of the film center. At steady state, the AV stays around 
the edge of PI without being fully engulfed. However, 
for (b), where the center of PI is (3/4) v2£ to the lower 
left of the center, the AV is completely absorbed into the 
pin. In Fig. 5 (c), the pins are bigger (d p = 2.5£) and 
their positions are as in (a). The AV is still not fully 
engulfed into PI, resulting in incomplete (hence weaker) 
pinning. We conclude that optimal pinning is sensitive 
to the precise location of the pins, and this sensitivity is 
not removed by using a larger pin. 

The dipole field strongly controls the positions of the 
V and the AV of the pair. Since the film is soaked in 
the applied magnetic field (Fig. 2), the V-AV pair is 
simultaneously under the influence of the dipole field and 
the pinning force. Thus, there should exist an optimal 
configuration of pin locations to obtain the best result 
for pinning. This optimum pin positioning is shown by 



Fig. 5 (b) (for the values of the parameters chosen) . 

V. MAGNETIC RELAXATION IN A FILM WITH 

PINS 

Due to the need for a larger computational domain, ex- 
cessive computational cost prohibits a complete simula- 
tion of magnetic penetration followed by magnetic relax- 
ation for larger V-AV separation. Since penetration and 
relaxation are two distinct processes from each other, we 
focus on relaxation to see if a V-AV pair can be pinned 
at the given separation. To that end, we take a single 
V-AV state from which the system can start to relax. 

Pinning of a V-AV pair is more difficult when the 
dipole field is turned off. The V-AV pair escapes from 
the pins and annihilates itself, even for the pin distance of 
(17/4)^2^ of Fig. 5 (b), which is the optimal (and maxi- 
mum) pinning distance we have obtained in the presence 
of the dipole field of <E> + = 1.90. Thus, based on obser- 
vations in previous figures, we have placed two pins at 
diagonal positions (with varying distances), and, follow- 
ing Ref. 18, an artificial V-AV pair. (See also 25-27 .) 

In several numerical experiments with random posi- 
tions of the V-AV pair as initial conditions, the steady 
states obtained are always the same as if the Meissner- 
state initial condition is used. This path-independence is 
expected for a thermodynamic equilibrium state. Thus, 
we have been able to create initial conditions in which a 
V-AV pair resides inside pins with arbitrary separation 
along the diagonal. By increasing the distance between 
the pins (between the V and the AV) we can reduce the 
attractive force between them. 15 

In Fig. 6 the pin distance is about 12\/2£ in a computa- 
tional domain of size 32£ x 32£. The system has reached 
a steady-state without showing any further movement 
of the V-AV pair from the pins. The pair initially at- 
tempted to escape the pins with a tendency to annihilate 
itself, but it eventually failed to emerge from the pins. 

Thus, we have achieved our goals (2) and (3) of the 
introduction to obtain a vortex state in which a V-AV 
pair is trapped in the pins, and is unable to overcome 
the pinning force and annihilate itself. 

Core deformation of the pair shown in Fig. 6 is char- 
acteristic of TDGL analysis (not observable in London 
analysis) such as Ref. 14. So the pin distance 12-\/2£ 
must be close to the critical distance that pins the pair 
at the pinning spots. (We observed annihilation of the 
pair at the pin distance of 8-\/2£, but did not systemati- 
cally search for the precise critical distance because each 
computation takes excessive computational cost.) 

If the V-AV pair remain pinned for a macroscopically 
long time after the dipole field has been removed, then 
taken together they can represent a 1 for nonvolatile 
data-storage, whereas the absence of a V-AV pair at the 
pins can represent a 0. In this case the density of vortices 
will be (1 V-AV pair)/(32£) 2 = 4.7GB /cm 2 , assuming 
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(a)Cooperpairdensity (b) Phase of older parameter 

FIG. 6. (Color online) The pins still hold the V-AV pair after a 
long time, in a 32£ X 32£ sample. We conclude that this computer 
modeling has reached steady-state and that pinning of the V-AV 
pair has been achieved. 

£ = 16 A for high temperature cuprate superconduc- 
tors. 28 One concern is the possibility that thermal fluc- 
tuations can overcome the energy barrier and de-pin the 
V-AV pair, which would cause our memory idea to fail. 
Within the Ginzburg-Landau theory we have estimated 
that the energy barrier of a single V-AV pair is indeed 
slightly lower than the energy level of thermal fluctua- 
tion (k B T) for YBa 2 Cu 3 7 ^s (YBCO) at around 77 K 
assuming d = 0.1£. However, the reverse can be true 
at somewhat lower temperatures (say, 0.5T C ). Alterna- 
tively, the barrier energy can be made larger by using a 
larger d, allowing our idea to become valid at, say 77 K, 
but the present theory is limited to d -C £, to allow us 
to perform a two dimensional calculation. Generalization 
to a three dimensional calculation is conceptually simple, 
but expensive in computational resources. 

So far, only the write part of this new storage device 
has been investigated, but we envision that reading can 
be achieved with a scanning magnetic force probe, a scan- 
ning Hall probe, or a scanning tunneling probe. Also we 
remark that by tilting the magnetic dipolc it should be 
easier to create a V-AV pair with the needed separation 
for pinning after the dipole field is removed. 

Our method to obtain precisely controlled vortex con- 
figurations may provide for wider, and more versatile ap- 
plications, e.g., Ref. 29 used a scanning SQUID micro- 
scope to create and move a single V-AV pair. In Ref. 29 
the field coil of the microscope was of 21/zm across, and 
hence applied a force over hundreds of square microns. 
In contrast, our method provides a controlled domain of 
submicron scale (only tens of coherence lengths) even at 
this early research stage. 

VI. SUMMARY AND CONCLUSION 

A numerical study of the magnetization process in a 
thin type-II superconducting film has been performed. 
The film is subject to a non-uniform magnetic field cre- 
ated by a dipole above the film. The dipole axis passes 
through the center of the film in the normal direction 



to the film plane. For a film with constant thickness 
and with no pin, it has been found that the film carries 
two V-AV pairs at steady state in a range of the im- 
posed magnetic flux <E> + (2.10 < $ + < 3.0), in units of 
flux quantum. The 3-pair-to-2-pair transitions, in which 
a pair of V-AV annihilates, have been observed in this 
range of <1> + . For <i> + < 2.10 no V-AV pair was created in 
the film. Two antidots (holes) have been used to create 
a single V-AV pair in the film, and a V-AV pair was cre- 
ated for lower magnetic fluxes to a minimum $ + = 1.3. 
It was observed that the balance between the dipole field 
and pinning force determines the optimal locations of 
the V and AV. A simulation of the magnetization pro- 
cess in a sample of size 32£ x 32£ suggests that we are 
likely to achieve pinning of the V-AV pair after removal 
of the dipole field with the V-AV separation indicated 
by this sample size. Using this sample as a template, 
the maximum density of pinned V-AV pairs achievable is 
calculated to be about 4.7GB /cm 2 . 
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